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Introduction 


• Research background: PWGL 

• Musical programming environment by Laurson, 
Kuuskankare, Norilo, Sprotte 

• High level visual interface to LISP programming 

• Synthesizer component in C++ written by the author: 
PWGLSynth 

• Kronos began as a bunch of aspirations for PWGLSynth 2 

• Generic computation engine 

• High level abstraction 

• Great performance 

• Since then, Kronos has morphed into a standalone 
compiler/language 

• Doctoral study project since 2010 
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The Current State of DSP Programming 


• The industry standard for DSP is C 

• Plugins and DSP chips tend to have C toolchains 

• AU, VST, LADSPA 

• Motorola 56k 

• etc.. 

• C is relatively hostile to casual programmers 

• Casual programmers make a lot of musical applications! 

• Getting audio out of C is very difficult for learners 

• High performance programs are low level 

• Many powerful abstractions have performance penalties 

• Tedious to write for professionals 
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Why Yet Another Programming Language? 

OR what if it would be possible to... 


• easily learn an audio language 

• write abstract, reusable code that runs fast 

• design all your algorithms down to the arithmetic primitive 

• have a single filter for any combination of single or double 
precision, real or complex, mono or multichannel... 


Many of us here are working on a subset of these problems. 
The final solution is not yet here. 
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A Language Specification 
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A Language Specification 

Syntax 


• Simple syntax 

• Familiar function notation 

SomeFunction(paraml param2) 

• Infix functions for arithmetics 

a + b * 3 / Sqrt(c) 

• Algebraic data structure yields pairs, lists and trees 

list = (abed) 

• Tie-in allows for partial decomposition too 

(first-element other-elements) = list 
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A Language Specification 

Classification 


• Functional Programming 
• No state 
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• Powerful abstraction 
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A Language Specification 

Classification 


• Functional Programming 

• No state 

• No variables 

• Similar to audio signal routing 

• Powerful abstraction 

• Reactive Paradigm 

• Action is followed by reaction 

• Reactive graphs are used to optimize signal rates 

• Implicit inferral of control and audio signals 
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A Language Specification 

An Example 


Listing 1: Fold, a higher order function for reducing lists with 
example replies. 

Fold(folding-function x) 

{ 

Fold = x 

} 

Fold(folding-function x xs) 

{ 

Fold = Eval(folding-function x Fold(folding-function xs)) 

} 

/* Add several numbers */ 

Fold(Add 1234) =>10 
/* Multiply several numbers */ 

Fold(Mul 5 6 10) => 300 
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• Kronos compiles programs on the fly to native x86 
• Some enhancements to SoftWire, a LGPL dynamic 
assembler written by Nicolas Capens 
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A Just-in-Time Compiler 


• Kronos compiles programs on the fly to native x86 

• Some enhancements to SoftWire, a LGPL dynamic 
assembler written by Nicolas Capens 

• Programs are configured on the fly for the present I/O 
configuration 

• Good runtime performance using SSE4.2 

• Comparable and often superior to an optimizing C-compiler 
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• Type Determinism 

• The result type of an expression can only depend on 
argument type 

• Complete dataflow analysis of the entire program 

• Result type can’t depend on data! 

• No dynamic containers! 

• Minimal branching! 

• No good for writing a word processor 
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Type Determinism 

...or The Catch 


• ...too good to be true? 

• To achieve performance, a rigorous constraint; 

• Type Determinism 

• The result type of an expression can only depend on 
argument type 

• Complete dataflow analysis of the entire program 

• Result type can’t depend on data! 

• No dynamic containers! 

• Minimal branching! 

• No good for writing a word processor 

• Designed for DSP inner loops 
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Signal Generation 


• Delay is a first-class unary operator 

• Recursion permitted via delays 

• Recursion can be turned into an osc by clocking a section 
of the loop 

• IO:Audio-Gen(sig) provides signal updates at the audio rate 
Listing 2: A Simple Phasor Oscillator 

Phasor(freq) 

{ 

next-phase = IO:Audio-Gen(z-1('0 phase + freq)) 
phase = next-phase - Truncate(next-phase) 

Phasor = phase + phase - 1 
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Fast -cos(7rx) for xe[-1,1] 


• Source is the math definition of series approximation! 
Listing 3: Taylor-series Cosine 


Sine-Coef(n) 

{ 

Sine-Coef = Crt:pow(Pi n * #2 - #1) * Crt:pow(#-l n + #1) / Factorial(n * #2 - #1) 

} 


Fast-Cos(x order) 

{ 

Use Algorithm 
xm = Abs(x) 
xp = xm - 0.5 

coefs = Map(Sine-Coef Count-To(order)) 

Fast-Cos = xp * Horner-Polynomial(xp * xp Reverse(coefs)) 
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Sinusoid Oscillator 


• Combining the - cos(tcx ) mapper and the phasor, a sine 
oscillator is created 

• With higher order functions, oscillator banks can be 
constructed from a list of frequencies! 

Listing 4: Sinusoid oscillation by mapping the phasor 

FSin(freq) 

{ 

FSin = Fast-Cos(Phasor(freq) #8) 

} 

/* Example of using higher order functions */ 

Osc-Bank = Reduce(Add Map(FSin freql freq2 freq3 freq4)) 
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SlN-OSC APPLICATIONS 


• Additive and FM Synthesis are easily constructed 


Listing 5: Sinusoid Synthesis 

Suppress-Alias(f0 amp) {Suppress-Alias = (fO <0.4) & amp} 

Additive(fO num-harmonics harmonic-coef harmonic-spread) 

{ 

Use Algorithm 

freqs = Map(Curry(Mul fO) Expand(num-harmonics Curry(Add harmonic-spread) 1)) 
amps = Expand(num-harmonics Curry(Mul harmonic-coef) 1) 

oscs = Zip-With(Mul(Map(FSin freqs) Zip-With(Suppress-Alias freqs amps))) 
Additive = Reduce(Add oscs) 

} 

FM(fO ratio mod feedback) 

{ 

modulator = FSin(fO * ratio) * mod 

carrier = FSin((l + modulator + feedback * z-l('0 carrier)) * fO) 

FM = carrier 

} 
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Interfacing with Control 


• Other clocks besides audio generators can be used 

• Resonator coefficient computations are clocked from OSC 

• Only recomputes coefs when OSC signal arrives! 

Listing 6: OSC-controlled Resonator Bank 

Reson(xO freq reson) 

{ 

xl = z-1 (' 0 xO) x2 = z-1('0 xl) yl = z-l('0 yO) y2 = z-l('0 yl) 
r = Crt:pow(reson 0.125) 

yO = xO - x2 + yl * 2 * r * Crt:cos(freq) - y2 * r * r 
Reson = yO * 0.5 * (1 - r * r) 


Reson-Bank() 

{ 

Use Algorithm 

params = ((IO:OSC-Input("cutoffl" Float) IO:OSC-Input("resonl” Float)) 
(IO:OSC-Input("cutoff2" Float) 10:OSC-Input("reson2" Float)) 
(IO:OSC-Input("cutoff3" Float) 10:OSC-Input("reson3" Float))) 
Reson-Bank = Reduce (Add Map (Curry (Reson NoiseO) params)) 
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SCHROEDER REVERBERATION 


• Classic Schroeder reverberation can be concisely 
expressed 

Listing 7: Classical Schroeder Reverb 

Feedback-for-RT60(rt60 delay) 

{ 

Feedback-for-RT60 = Crt:pow(#0.001 delay / rt60) 

} 


Basic(sig rt60) 

{ 

Use Algorithm 

allpass-params = ((0.7 #221) (0.7 #75)) 

delay-times = (#1310 #1636 #1813 #1927) 

feedbacks = Map(Curry(Feedback-for-RT60 rt60) delay-times) 

comb-section = Reduce(Add Zip-With(Curry(Delay sig) feedbacks delay-times)) 
Basic = Cascade(Allpass-Comb comb-section allpass-params) 
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FDN Reverberation 


• Some highlights from a FDN reverberator 

Listing 8: Snippets of a 16th order FDN reverberator 


/* Orthogonal matrix multiply - Householder Algorithm */ 

Feedback-Mtx(input) 

{ 

Use Algorithm 
Feedback-Mtx = input 
(even odd) = Split(input) 

even-mtx = Recur(even) odd-mtx = Recur(odd) 

Feedback-Mtx = Append(Zip-With(Add even-mtx odd-mtx) Zip-With(Sub even-mtx 


odd-mtx)) 


/* 16-channel feedback signal recursively passed through a unit delay operator */ 
feedback-vector = z-l(' (0 00000000000000 0) 


Zip-With(Mul loss-coefs Zip-With(Filter:OnePole 
Feedback-Mtx(delay-vector) filter-coefs))) 
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Thank you! 


...questions? 
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